Fracture and second-order phase transitions. 



O 

o 
o 

(N 

< 

00 



I 

O 

o 



> 

00 

O 
O 
O 

> 



I 

O 

o 



X 
S3 



Y. Moreno^y, J. B. Gomez^, A. F. Pacheco^ 
^ Departamento de Fisica Teorica, Universidad de Zaragoza, 50009 Zaragoza, Spain. 
^ Departamento de Ciencias de la Tierra, Universidad de Zaragoza, 50009 Zaragoza, Spain. 

(February 1, 2008) 

Using the global fiber bundle model as a tractable scheme of progressive fracture in heterogeneous 
materials, we define the branching ratio in avalanches as a suitable order parameter to clarify the 
order of the phase transition occurring at the collapse of the system. The model is analyzed using 
a probabilistic approach suited to smooth fluctuations. The branching ratio shows a behavior 
analogous to the magnetization in known magnetic systems with 2nd-order phase transitions. We 
obtain a universal critical exponent ~ 0.5 independent of the probability distribution used to 
assign the strengths of individual fibers. 
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The interest in the fracture processes of heterogeneous 
media has increased in the last years In the lab, 

a disordered material subjected to an increasing external 
load can be studied by measuring the acoustic emissions 
before the global rupture. It has been shown that 
this intense precursory activity in the form of bursts of 
different microscopic sizes follows a well-defined power 
law. In despite of the many efforts and successes that 
have been recently achieved, the question of whether rup- 
ture exhibits the properties of a first-order or a second- 
order phase transition remains under discussion as well 
as what is the order parameter that indicates the type of 
transition. 

From the theory side, the understanding of fracture 
in heterogeneous materials has progressed due to the 
use of lattice models and large scale simulations |3| . In 
this field, it is important to use models able to describe 
the complexity of the rupture process, nevertheless, they 
should be simple enough as to permit analytical insights. 
To this class of models belong the well-known Fiber Bun- 
dle Models (FBM) widely used since their introduction 
more than forty years ago 0,^. In static FBM, a set of 
fibers (elements) is located on a supporting lattice and 
one assigns to its elements a random strength thresh- 
old sampled from a probability distribution. The set is 
loaded and fibers break when their loads exceed their 
threshold values. In the Equal Load Sharing (or global) 
FBM, which is the simplest scheme one can adopt to 
make the problem analytically tractable, one assumes 
that the load carried by failed elements is equally dis- 
tributed among the surviving elements of the system. 

In the present Letter, we explore the criticality of frac- 
ture in the global FBM using a novel probabilistic ap- 
proach devised to smooth fluctuations. The scaling rela- 
tions obtained and the behavior of the order parameter 
point out that the fracture of a fiber set with long range 
interaction undergoes a continuous phase transition. 

Let us first recall the basic ingredients of the global 
FBM and how one proceeds in numerical simulations. 



The system under consideration is a set of Nq elements 
located in a supporting lattice each one having at the 
initial state a zero load and a fixed strength threshold 
value sampled randomly from a probability distribution 
P{<j). The system is then subjected to an external force 
F such that each element increases its load in the same 
amount, a. This individual stress (or load), a, will act 
as the control parameter. The process of driving is done 
quasi-statically; i.e., the external force is increased at a 
sufficiently slow rate as to produce a single breaking event 
when the stress on the weakest element equals its thresh- 
old value. Then, the increase of F stops and the load of 
the broken element is equally transferred. This implies 
that the load on any element is given by ti = F/na{F), 
where Us is the number of surviving elements for a given 
F . The rupture of an element may induce secondary fail- 
ures which in turn may trigger more failures and so on. 
This process of induced failure at constant external load, 
termed an avalanche, stops when all surviving elements 
carry a load lower than their thresholds. The system is 
then loaded again and the process is repeated until the 
final catastrophic avalanche provokes the total rupture 
of the material, which occurs at a critical load Uc that 
depends on the probability distribution from where the 
individual strengths were drawn, as well as on the system 
size. The FBM have been recently used in self-organized 
criticality (SOC), a theoretical framework widely used 
for the study of avalanche phenomena in disordered sys- 
tems. It has been shown using these models that systems 
with plastic behavior can reach a SOC state just before 
the global rupture [|j . A second case of self-organization 
with power law distributions in several quantities cor- 
responds to the situation in which the fracture process 
coexists with a healing process ||T^ . 

In numerical simulations, the cycle of complete break- 
down of the model is performed many times in order to 
average out the effect of fiuctuations and obtain mean 
values. As we are interested in studying the behavior of 
the system as the critical point, or point of final collapse, 
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is approached, it is of most importance to find a sim- 
ple method able to capture the evolution of the system 
avoiding as much as possible the fluctuations appearing 
in numerical simulations. To introduce our probabilistic 
strategy, let there be a large set of iVo elements. Suppose 
that each element carries a given load cr, which is zero at 
the initial state. The strength of each element is drawn 
from a probability distribution P{<y). Different probabil- 
ity distributions can be considered. In materials science 
the Weibull distribution is widely used. 



P{a) = 1 - e 



V crn ' 



(1) 



p being the so-called Weibull index, which controls the 
degree of disorder in the system (the bigger the Weibull 
index, the smaller the disorder), and (Tq is a load of ref- 
erence. In the following we will assume ctq = 1, and 
therefore the loads are dimensionless. At this point, it 
is worth noting that the results and the formulas de- 
rived in the following hold for a wide class of probabil- 
ity distributions. We use here the Weibull distribution 
for definiteness but results have also been obtained for 
other distributions. Equation ([l]) represents the proba- 
bility that an element fails under the individual load a. 
Now, consider the case in which an element drawn from 
Eq. (Q) supports a load ci but breaks under a new load 
a2- The probability that this happens is given by 



p(o-i,cr2) = 



P{a2) - Pjai) 
l-P(ai) 



(2) 



So, the probability q{(7i,(72) that an element that has 
survived to the load ai also survives to the load (T2 will 
be given by g(cri, (T2) = 1 — 0-2) = e~^'^2-'^i). 

To mimic the quasi-static increase in load on the sys- 
tem we impose the condition that under an external force 
F, the next breaking event consists of one single fail- 
ure. Let suppose that after the latest avalanche, there 
are Nk surviving elements each one bearing a load Cfc. 
The new individual load ai needed to provoke the fail- 
ure of just one more element is given by the solution of 
Nk-l = Nk-q{cTk,cTi). Thus, 
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where in Eq. (^ N-k = A'o and CTfe = at the initial state. 
Elevating the external force up to the N]^ ■ ui level, statis- 
tically speaking, one element breaks. As we are dealing 
with an equal load sharing set, the choice of the broken 
element is irrelevant because all the elements are equiva- 
lent. Once the first element fails, the redistribution of its 
stress takes place. This may induce other failures until 
the end of the avalanche. 

Now, how many elements will survive to the situation 
in which n\ elements with load a\ fail in an avalanche 



step? The new load on the intact N\ — n\ = N2 ele- 
ments is (72 = ^^4^^ ■ So, the number A^3 of elements that 
survive to the new load can be expressed as 

N-i^N2- q{(Ji, — • CTi) = iV2 • g(cri,CT2). (4) 

1\2 

As a consequence of applying Eq. (^, A^2 " N3 elements 
break and the new total number of intact fibers will sup- 
port a bigger load 173. The avalanche may continue and 
Eq. (^ is applied again for the set of surviving el- 
ements. The iterative process will stop when no new 
element fails, which occurs when the right-hand side is 
equal to the left-hand side in Eq. (^). The general form 
of Eq. (H) is 



(5) 



with the conservation condition for the total load in the 
system during an avalanche 



and the condition 



(6) 



(7) 



which determines the end of the avalanche. The dynam- 
ics of the system is determined by Eq. (^), (||), (^. The 
size of an avalanche is given by the number of elements 
that break between two successive steps of external load- 
ing. The critical load, defined as the load needed to pro- 
voke the total collapse of the system, is equal to the load 
on the intact elements just before the final catastrophic 
avalanche. Note that in this probabilistic approach, in 
contrast to Monte Carlo simulations, we only need to 
store the information concerning the loads on the intact 
elements, that is, the details of the threshold list is omit- 
ted. 

In the probabilistic strategy, we can proceed in two 
different ways in order to determine when an avalanche 
ends, to which we will refer as the continuous and the dis- 
crete cases. For the continuous case, the number Nj^i of 
surviving elements is considered a real number. Strictly 
speaking, this means that the condition (Q) is never ful- 
filled. So, the condition (|^) is replaced in numerical cal- 
culations by using a factor <C 1 that determines the 
end of an avalanche, i.e., if A'^ — Nj^i < v the avalanche 
stops; otherwise it continues. In the discrete case, Aj+i 
is considered to be a whole number, so that after each 
iteration of Eq. (||), A^+i has to be rounded up. This 
is done comparing the remainder of A'j+i, A, with a ran- 
dom number a uniformly distributed in the interval [0, 1[. 
Thus, if a > A, A^+i is equal to its whole part, other- 
wise, ATj+i is equal to its whole part plus one. Next, 
we check whether the condition (^ is satisfied for the 
rounded value of N j+\ or if a new iteration of Eq. (|^) 
has to be performed. The continuous approach has the 
advantage that the fluctuations are ruled out, whereas 
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for the discrete case the results are similar to those ob- 
tained by Monte Carlo simulations where it is necessary 
to average over many realizations in order to get accurate 
mean values. Remember that in this model the central 
limit theorem applies 

In Fig. d we have depicted the fraction of broken el- 
ements versus cr, for the continuous case and for four 
individual Monte Carlo simulations with a WeibuU in- 
dex p — 2 and A^'o = 5000. No averaging has been done 
because our aim is just to illustrate the scatter of the re- 
sults. As can be seen, the continuous probabilistic model 
gives a smooth curve, and provides an accurate value 
for the critical load CTc, which analytically is given by 
CTc = {pe)-'^/P [D in the hmit of infinite TVq. 

Now, we proceed to explore the behavior of some quan- 
tities as the critical point is reached. The results shown 
below have been obtained using the continuous approach 
(p = 2). In Fig. H we show the interesting scahng rela- 
tion for the average avalanche size. It turns out that the 
avalanche size near to the critical point diverges with an 



exponent 7 



as s 



A similar behavior, 



through a mapping of a fuse network model to the global 
fiber bundle model used here, has been recently reported 
. We have also obtained the same scaling function for 
the derivative of the number of broken fibers with respect 
to the load on the system. The rate dN/da diverges as 
(ctc — (t)~ 2 ^ thus qualifying a critical mean field behavior 
as was already shown in Ref. [Q. In Ref. a simi- 
lar scaling behavior is addressed for the derivative of the 
strain carried by the fibers with respect to the driving 
force. 

Another way to shed light on the critical behavior of 
this type of system is to define a branching ratio ( for 
each avalanche. This magnitude represents the proba- 
bility to trigger future breaking events given an initial 
individual failure [ p^[p^ ] , and is related to the number of 
broken fibers by 



c 



< z > -1 
< z> 



(8) 



The above relation can be obtained by thinking of the 
evolution of fracture as a kind of branching process [ p^ . 
In this process, each node gives rise to a number n of 
new branches in the next time step. The average number 
< n > of new branches is called the branching ratio. Let 
us denote by nt the number of branches at a given step t 
of the branching process, and by tmax the total number 



of time steps before it dies. Then, C = 1 



As, 



no = 1, C = 1 ~ where ritot is the total number of 
nodes developed in the branching process. For a fracture 
process, Utot is equal to the average number of failure 
events. So, Eq. defines the branching ratio. We rep- 
resent by < z > the average number of elements that fail 
in one avalanche, which is a function of a and coincides 
with s. This analogy between fracture and branching 



processes has been previously used to study the critical- 
ity in the process of fragmentation of Hg drops . The 
branching ratio will then act as the order parameter. It 
takes the value 1 when the system is critical thereby rep- 
resenting a measure of the distance of the system from 
the critical state Q. We have numerically computed C 
by means of the continuous method. The results obtained 
for a system of A^o = 50000 elements and several values 
of p have been plotted in Fig. |3[ It can be seen in this 
figure that 1 — C approaches zero, in all cases, in a con- 
tinuous fashion as the critical load is reached. Near the 
critical point, the relation 1 — ( ^ {(Tc — cr)^, with f3 — ^ 
applies. This is the exponent of the order parameter 
that we should expect from a mean-field approach. Note 
the similarity of the figure with those obtained for the 
magnetization in known magnetic systems with second- 
order phase transitions. In the figure, the values of C are 
collected for all the avalanches except for the last, that 
which provokes the collapse of the system. The result 
that at (Tc C ~* 1; is consistent with the previous result 
that the avalanche size diverges at the critical point. On 
the other hand, the branching ratio does not depend on 
the size of the system for large systems, in contrast with 
previous results in other fracturing systems p^ . 

The suggestion of Ref. |^ is that fracture can be seen 
as a first-order phase transition close to a spinodal-like 
instability |16|. There, by simulating models of electric 



breakdown and fracture, the authors present numerical 
and theoretical evidence of several scaling relations and 
of a discrete jump in some macroscopic properties. Here, 
we have obtained the same scaling relation for the rate of 
fiber failures as the critical point is reached, and for the 
avalanche sizes, which also diverge at that point. Our nu- 
merical results also fit the mean field result 7 = 5- It is 
true that the fraction of unbroken fibers just before the 
global rupture has a discontinuity; but from our point 
of view that is not enough as to set the conclusion that 
fracture can be describe as a first-order phase transition, 
since the concepts related to spinodal nucleation are not 
sufficiently well established in driven disordered systems. 

Our alternative point of view has been to consider 
the branching ratio as an appropriate order parameter. 
According to the results obtained, the branching ratio 
goes continuously from zero to one. Note, additionally, 
that what changes discontinuously at cTc is the rate of 
change of ( rather than ( itself, which is in the essence 
of a continuous phase change. Therefore, the behavior of 
the branching ratio implies that the system undergoes a 
second-order phase transition as claimed in other analysis 
of fracture models . Our results suggest that fracture 
in heterogeneous systems with long range interaction can 
be described as a phase transition of the second-order 
type. Besides, it is important to bear in mind that in 
fiber models with local interactions, the order parameter 

has a discontinuous jump typical of first order phase 
transitions. Finally, it is worth recalling that fracture of 
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real materials is a process based on elasticity, and elas- 
ticity is a long-range phenomenon. In this respect, the 
global load-sharing model we have explored here could 
be a better analogy to real fracture than the local one. 
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FIG. 1. Fraction of broken elements for the equal load shar- 
ing model. The line corresponds to the results obtained with 
the continuous approach and gray dots correspond to four 
Monte Carlo realizations. 




FIG. 2. Scaling of the mean avalanche size, s, as the crit- 
ical point Oc is approached. The results correspond to the 
continuous probabilistic approach for a system of A^o = 50000 
elements and p = 2. The straight line with a slope — ^ has 
been drawn for comparison. 
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FIG. 3. Evolution of the branching ratio as the critical 
point is approached in the continuous probabilistic method 

{No = 50000). Note that at the critical point the branching 
ratio reaches the unity. The critical exponent is /9 = i which 
coincides, as it should, with the value of the exponent 7. 
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